Critical wetting of a class of nonequilibrium interfaces: a mean-field picture 
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A self-consistent mean-field method is used to study critical wetting transitions under nonequilib- 
rium conditions by analyzing Kardar-Parisi-Zhang (KPZ) interfaces in the presence of a bounding 
substrate. In the case of positive KPZ nonlinearity a single (Gaussian) regime is found. On the 
contrary, interfaces corresponding to negative nonlinearities lead to three different regimes of criti- 
cal behavior for the surface order-parameter: (i) a trivial Gaussian regime, (ii) a weak-fluctuation 
regime with a trivially located critical point and nontrivial exponents, and (iii) a highly non-trivial 
strong-fluctuation regime, for which we provide a full solution by finding the zeros of parabolic- 
cylinder functions. These analytical results are also verified by solving numerically the self-consistent 
equation in each case. Analogies with and differences from equilibrium critical wetting as well as 
nonequilibrium complete wetting are also discussed. 

PACS numbers: 02.50.Ey,05.50.+q,64.60.-i 



I. INTRODUCTION 

Wetting transitions can be considered as an interface unbinding transition where the mean interfacial separation 
from a substrate plays the role of an order-parameter. In the wet state the interface drifts away from the substrate, 
while in the nonwet state the interface is bound to it, with many contact points between them. In equilibrium, such 
phenomena can be described by Langevin-type equations of the form [l[ 



dh(x,t) nr7j , SV(h) . , 

with V(h) standing for the energy between the substrate and the interface, r\ being a Gaussian white noise, and D 
a diffusion coefficient. Alternatively, the number of contact points between the substrate and the interface can also 
be regarded as an order-parameter that equals zero when (h) = oo and is nonzero otherwise. This quantity is closely 
related to the surface order-parameters studied in the framework of lattice systems in a semi-infinite geometry, for 
which the interfacial displacement models like Eq. (JTJ are expected to be useful effective descriptions [2] . A familiar 
example is provided by Ising ferromagnets where the surface order-parameter is the average magnetization at the 
substrate. The variable (n) = (e~ h ) constitutes an adequate mathematical representation of such an order-parameter 
that exhibits singular behavior: it is a positive quantity that vanishes as (n) ~ \a — a c \^ close to criticality, with a 
denoting a convenient control parameter, and evolves in time as (n) ~ t~ e right at the transition. 

An extension to nonequilibrium situations was considered by adding the nonlinear term A(V/i) 2 to Eq. ([T]), which 
represents preferential growth along the local normal to the surface and covers the realm of Kardar-Parisi-Zhang (KPZ) 
interfacial phenomena [1, 0]. This is expected to capture the physics of wetting transitions under nonequilibrium 
circumstances in the simplest possible form. Interestingly, in this case a simple Langevin equation for n does exist if 
short-ranged forces between the substrate and the interface are assumed. One typically writes [l|, 

V(h) = [ dx \ah(x) + be- ph W] (2) 



where a and b are phcnomcnological parameters. The change of variables h = sgn(A) Inn leads to an equation of the 
form (see [j| for details), 



M*, t) = DV 2 n + an + bn i+ P + / ~ / 3 n 
dt 

to be interpreted in the Stratonovich sense. Here, & is a temperature dependent-parameter, a plays the role of a 
chemical potential difference and is the control parameter, and (r/(x,t)) = 0, (t](x, t)rj(x' , t')) — a 2 5(x — x')S(t — t'). 
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p carries the opposite sign of A which is the ultimate responsible for the physical behavior. For p > 0, the above 
equation represents the so called multiplicative noise 1 (MN1) universality class [H H[ . For p < 0, n is merely an 
auxiliary variable that diverges at the transition and it is 1/n which has to be studied instead. In this case the 
associated set of exponents define the multiplicative noise 2 (MN2) universality class @,0,[1]- A detailed discussion 
of the differences between the MN1 and MN2 universality classes can be found in Q . In the following discussion the 
surface order-parameter will be denoted nop- It should be emphasized that nop = n for MN1 and nop = 1/n for 
MN2. Despite the simple appearance of Eq. ([3]), we nonetheless caution that the intricacies and subtleties of the KPZ 
equation lie in the interplay between the diffusion and the multiplicative Gaussian white noise terms. 

The MN1 and MN2 universality classes are examples of nonequilibrium, complete wetting transitions, in which a has 
to be fine tuned. A second type of wetting transitions may occur if as the temperature is increased b becomes bigger 
and eventually changes its sign while the system is kept at coexistence, a — a c . This is denoted critical wetting and 
amounts to taking b < and adding a higher-order term cn 1+2p in order to obtain a finite solution for n. Equivalcntly, 
the potential to be considered in the interfacial representation is V(h) = ah + be~ ph + ce~ 2ph , where the last term 
ensures stability. More precisely, the model system is defined by the Langevin equation 



^' ' = DV 2 n + an + bn 1+p + cn 1+2p + nn{x, t), (4) 

in which now b is the new control parameter and a is set to the critical value a c found for the complete wetting 
transition. For sufficiently low values of b < the interface remains pinned and the density of locally pinned sites 
at the wall is high. As the transition is approached (increasing 6), the stationary density of pinned segments goes to 
zero in a continuous manner as (nop(b,t — oo)) ~ \b — b c \". Above b c , the interface depins and therefore the mean 
separation (h) diverges and {nop} vanishes. In this latter case the density of pinned segments scales with time as 
(n p{b > b c , t)) ~ t~ 6 ', with the exponent 9 adopting different values for MN1 (p > 0) and MN2 (p < 0). 

In this paper we study within the mean-field approximation the critical wetting transition associated with Eq. 
(j4|). In the absence of exact solutions, mean-field approaches are useful not only in enabling analytic calculations 
to be performed, but also because they provide insight into the physical behavior at high system dimensionalities 
which would be otherwise unattainable from computer simulations alone. Since, in the present case, it is known that 
mean-field theory could be valid for dimensions as low as d — 2 (d — 3 bulk dimensions) at least in some regime, the 
results presented here may be relevant for realistic three-dimensional systems [5| . 

This paper is organized as follows. In Sec. II we describe the mean-field approach and provide details of the 
calculation for both the MN1 and MN2 cases. As will be proved, three different scaling regimes have to be distinguished 
for MN1, and only one for MN2. Results for higher moments of (nop) and for (h) are also included. Section III 
contains a summary and a discussion of our findings. 



II. MEAN-FIELD APPROACH 



Sound mean-field approximations to multiplicative-noise equations like (U]) require that the effects of both the noise 
and of the spatially varying order-parameter are taken into account to some extent. To this effect, the following 
procedure can be used [9(: the Laplacian is discretized as l/2dj^ .(rij ~ n i)i where rij = n(xi, t) and the sum is over 
the nearest- neighbors of i. Afterward, the value of the nearest-neighbor is substituted by the average field (n) to obtain 
a closed Fokker-Planck equation for P(n,t, (n)). The steady-state solution is then found from the self-consistency 
requirement [l(| 



<„> . f" P '"' <"» (5) 
Jo p (n,{n)) 

In what follows we particularize Eq. © to the MN1 and MN2 cases [see Eq. Q]. It will be shown that for MN1 
there is a sequence of scaling regimes depending on the relative importance of the noise strength as compared to the 
spatial coupling D and the nonlincarity exponent p: (i) a pure mean-field regime where the noise can be completely 
disregarded, (ii) a weak-noise regime where the noise strength enters the expression of the exponents, but without 
shifting the wetting temperature, and (iii) a strong noise regime where both the wetting temperature and the exponents 
are noise dependent. For MN2 the situation is far less rich, exhibiting a single mean-field-like scaling regime. 

In the following analysis the mean- value theorem for infinite integrals [ll| will be used repeatedly to determine the 
asymptotic behavior in (n) of the various integrals. According to this theorem, under quite general integrability and 
boundedness conditions. 



3 



dx f{x)g{x) = [i j dx g(x), (6) 

J a 

where fi is some value between the lower and upper bounds of f(x) [ll|]. Likewise, the combination 2D /a 2 will be 
denoted by v to simplify the notation. 

A. The case MN1 

For the MN1 case (p > 0), the associated stationary probability density can be readily obtained from the associated 
associated Fokker-Planck equation [l2j and reads 

P -(B)K ^. <Bp {-£(i B , + !»* + 2!)}, ( 7) 

where m = (n), and 7 = — 2 (a + D)/a 2 . After defining 



1 (ml = / an n' exp < =— rr ~— n ^ 5 — >, (8j 

7 [ cHp a A p a A n ) 

and substituting n = vmx, I{m) = (vm) 1+1 J (m) , with 

J = /°° dx ^e-'Vl'-cV) 8 ^-!/^ b > = ^ d = -^- v 2p . (9) 
Jo c p a z p 

The self-consistency equation ([5]) can now be recast in the simpler form 

v 1 + 7 , d m J(m) 



m m J(m) 



(10) 



Since the mean-field, self-consistent calculation for complete wetting (MN1) [T(| yields a c = a 2 /2 (j c + 1 = v) and 
J > does not diverge, the condition (|10p simplifies to d m J = 0. 

We now consider an intermediate point x\ > such that mx\ <C 1 and split J(m) as 

rx\ poo 

J = Jl + J2 = / dx x -l-» e -b'('nxr-c'( m xf» e ~l/x + / dx x -l-v e -b>{mxy-c'{mxy e -l/ X ^ ^ 
JO Jxi 

after which exp[— b'(mx) p — c'(mx) 2p ] in J\ is expanded to second order, whereupon 

Ji ~ Cl - b'c 2 m p - c 3 m 2p + 0(m 3p ), (12) 
where c\, C2, and C3 are constants. As for J 2 , since the derivative can enter the integral, 

' dx x- 1 - v e- b ' (mxY - c ' {mx)2P e- 1/x (-pb'm p - 1 x p - 2pd m 2 ^ 1 x 2p ) 

POO 

= e -VeM / X "i-^ e -f''("^) p -c'(™^) 2p / -pl' m v-^ x v - 2pc'm 2p - 1 x 2p Y (13) 

Al 

with £(m) G [xi, oo) being a by-product of the application of the mean- value theorem. Taking (mx) p = t, 



dm 



dm m 1 v J/ mxi )%> V ) m 
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Finally, 

^- = pb'czmP- 1 + 2pc 3 m 2p - 1 + cAm)m v - 1 = 0. (15) 
dm 

To proceed further requires identifying the term involving the lowest power of 777, which in turn requires studying how 
04(171) modifies m v ~~ 1 for small values of m. 

First, notice that the factor exp(— 1/£) is innocuous in such a limit. Second, we work out the low-m limit of the 
two integrals contained in €4(771), namely, 



/■OO 

C= / dtr* e- bWt2 , Cu+x= dtt-f +1 e- b,t - c ' t2 , (16) 

(mxi) p J (mxi)P 



to find after splitting, 



/ dt i 5 e- bW ' 2 = / dtt 5 e- b,t - c ' e + / dtt s e- b,t - c ' t2 . (17) 
After applying the mean-value theorem to the first integral a contribution m p(i+<*) i s 

obtained (the second integral 

contributes a constant), whence it ensues that the leading asymptotic behavior of Eq. (|15|) is unaffected by 04(777) 
and hence two cases must be distinguished. 

Case 1: p < v leads to an equation of the form (02,03, and C4 being positive factors) 

bc 2 + m p c 3 + m^ci = 0, (18) 

which implies m ~ (— b)^ as b ^ 0, with 

if (19) 

Cast 2: p > v. It is expedient to rewrite Eq. (fT^|) as 



>2 

dm 



p(rnxi) p \ 

t -v/ P -i e -b't-c't* ^ t + 2c /^ _ j dt t -u/ P -i e - b 't-c'e ^ t + 2c , t 2^ j 



dt 



(20) 



The first integral, which we call 05(b), does not depend on m but can vanish for particular values of b, while the for 
the second one we again apply the mean- value theorem to obtain the leading (lowest) powers mP~ v ', m 2p ~ v when 
to ^ 0. Last, 



bc 2 mP- v + c z m 2p - v + cb(6) = 0. (21) 

Note that the factors £2 and 63 are given in terms of exponentials of £1 £ [x±, oo),^2 £ [0, (mx\) p ] which result from 
the application of the mean-value theorem, and do not vanish as 777 — + 0. The next step is to find out what values 
of b make 05(b) vanish. The latter consists of two integrals that can be easily written in terms of parabolic-cylinder 
functions D fJ/ (x), using the recurrence relation D /1 (x) = xD^^i(x) + (1 — (i) Dp— 2(1) [Tit : 



05(b) = (2,/)* r(l - £) e^D f k/^Q • (22) 

Therefore, the critical point is determined by the zeros of the parabolic-cylinder functions. It turns out that D v / p 
has exactly one zero of order one on the interval at hand 13], < vjp < 1, and that this zero is negative, so we 
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thus conclude that in this regime the transition occurs at the finite value b c = y ^-^-x c , with x c < the only zero of 

D v i p (x), and is controlled by an exponent (3 = l/(p — v), where now (n) ~ (— b + b c )P . 
We next summarize the scaling regimes obtained for the three cases: 



1 

V 



(b c = 0) if p<£ 1 , 
a 



2D 



1 - (6 C = 0) if 4 <P< ™, 
cr a 

P 



(23) 



(b c <0) if P > 2D m 



P 2 

In order to check these results we have solved numerically the self-consistent equation ([5])- This requires evaluating 
numerically the involved integrals. Figure 1 illustrates the output of this calculation by showing estimates of (3 as 
a function of p and 2D /a 2 . Note the excellent agreement with the analytical results Eq. (f2"3")) . In the region around 
a 2 p/2D w 1, there are divergences and the integrals are difficult to evaluate numerically, generating large error 
bars. For ratios a 2 pi 2D larger than 1 the location of the critical point obtained by solving numerically is find to 
coincide with that given by the zeros of the corresponding parabolic-cylinder function, which we have also computed 
numerically. These results provide a complete verification of the previous analytical calculations. 

Let us now consider higher-order moments. We can write for k > 



m k = (n k ) 



h{m) 
I (m) 



(24) 



with 



Jfc(m) 



dn n 



7+fe-i 



exp 



2b 



a 2 p 



a 2 p 



- ■ " n 



(25) 



where 



dx x^+k- 1 e - b ' ( mx ) p - c ' ( mx ) 2p e ~ 1 / x 



(26) 



Splitting the integral into two parts as above , — j[ k ^ + J. 



r(fc) 



J« = jW + J. 



(fc) 

2 



dx x ^+ k - 1 e ~ b '( mx ) p - c '( mx ) 2p e -i 



+ 



dx x~ /+k ~ 1 e~ b '( mx y '~ c 'i mx ) 2 



(27) 



Proceeding similarly as in the above cases for Ji and J2, the lowest powers in m can be identified as 



(k) 



Ik) 



b'^mP - c { 3 k) m 2 P + 



(28) 



and using ([17] 



t(*0 
J 2 



„(fc) 
c\ 



b'c^m-*-* 



(29) 



Consequently, 



TOfe 



if k< + 
if k> + 



(30) 



This represents a strong form of multiscaling, very similar to the one reported for the moments in the self-consistent 
solution for MN1 [l|. 

Also, (h) = (— Inn) can be computed effortlessly by making use of 



In n = lim 



n a - 1 



(31) 
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which reduces the calculation of the average of Inn to a combination of moments of n. Use of this gives (h) ~ 

-H-b + b c ). 

Finally, we just mention that the more general situation bn 1+p + cn 1+q with q > p leads to the following simple 
substitutions: 

q^p if P<9<7- 
^1 if p < 2D < ^ ) 

/3 = <( — - P !T (32) 

°" 1 if „>2£ 

To locate the critical point one has now to proceed numerically to search for the zeros of c$(b) which can no longer 
be expressed in terms of parabolic-cylinder functions. 

B. The case MN2 

Consider again Eq. ([¥]) where for convenience we have introduced a change in the sign of p, which is now positive, 

dn(x, t) 



i)t 



DV 2 n + an + bn 1 ^ + en 1 - 2 ? + nn(x, t). (33) 



This is a non-order-parameter Langevin equation whose associated universality class can be studied by measuring the 
order-parameter The corresponding stationary probability density is 

PAn) - exp { - A ( »„-p + + EHt)\ (34 ) 

{ <J \p 2p n J J 

with m = (n) — > oo at the transition. Proceeding as in the MN1 case, we define 

7(m) = dn exp (--?-»-> - -^n" 2 * - ^-V (35) 

Jo V °" P °" P a z n J 

and the self-consistent equation is now given by 

v _ d m I(m) 



m I(m) 



(36) 



This last expression can be further simplified by making the change of variable x = m/n and setting as before 
a c = a 2 /2 [8| to obtain 

flUJ(m) = Qj (3?) 
J(m) 

where 

poo 

I{m) = m- u J{m) 1 J{m) = dx x ^ e -b'^/m.r-c'(x/ m ) 2 " e -u/ x ^ ^ 



Jo 

and v, and c' are defined as above. Next, an intermediate point x\ is considered such that << 1 an d J(m) is 
split as in equation (jTTj) . J(m) = Ji + </2- The exponential factor exp[— b'(x/m) p — c'(x/m) 2p \ in J\(m) is expanded, 
what leads to the following asymptotic behavior 

Ji(m) ~ ci - b'c 2 m p - c 3 m 2p . (39) 
Regarding J 2 , by virtue of the mean- value theorem 

= -b'pm-^MO - c'2pm' 2p -\m) (40) 

dm 
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where /<(£) stands for the function c , ie - b 'K*(" l )/ m ] p -c%M/™] 2p w ith &( m ) e [zi.oo) and c, = /~ da; x a *e- ra Note 
that a, does not depend on m. Lastly, after substituting the self-consistent equation reads 

4- (ci - b'c 2 m-P + J 2 (m)) - 0, (41) 
dm 

from which it is straightforward to extract the critical point b c = and write 

m ~ (-b)- l ' p . (42) 
This result has been verified by numerical integration of the self-consistent equation. 

A calculation of the higher moments of the distribution (|3"4"|) along the same lines as the previous section (but with 
m — > oo ) results in mj- = (n fe ) ~ m k , and therefore {nop) ~ (— fo) 1 ^ while, as in MN1, (/i) ~ ln(m) ~ — ln(— b). 
Results pertaining to (nop) are numerically verified in Fig. 2. 



III. DISCUSSION AND SUMMARY 



We have investigated within the mean-field approximation the surface order-parameter n — e~ h at nonequilibrium, 
critical wetting transitions of KPZ interfaces interacting with a wall. The model, as described by the Langevin 
equation covers both positive (MN2 p < 0) and negative (MN1 p > 0) KPZ nonlinearities. 

The more interesting case is that of negative KPZ nonlinearities, where three different scaling regimes can be 
distinguished. 

1. If p < D/a 2 a critical wetting transition exists characterized by the exponent (3 — l/p and the critical temper- 
ature b c = 0. We denote this the pure mean-field regime because these are the values that would have obtained 
had the noise and the Laplacian terms been neglected in the equation, i.e. in the crudest possible mean-field 
limit flU). 

2. If D/a 2 < p < 2D I a 2 the critical wetting temperature is still given by b c = but the critical exponent no longer 
depends solely on the potential details but also on the noise strength, (3 = l/(2D/a 2 —p). This is denoted the 
weak-noise regime. 

3. If p > 2D /a 2 the system enters a strong-noise regime where the critical temperature is shifted away from zero, 
b c = x c \J a 2 pc/2 with x c < the only zero of the parabolic-cylinder function D 2 D/a 2 p( x ): and /3 = \ / [p~2D / a 2 ) . 

This rich structure is expected to hold at least qualitatively beyond mean-field theory as it is known that Eq. ([3]) 
exhibits a strong-coupling regime for arbitrary high system dimensionalities. A similar scenario arises in the full 
solution of equilibrium critical wetting with long-ranged forces for which there are also three regimes for the behavior 
of (h), whose nature is similar to ours insofar as their origin can be traced back to the relevance of fluctuations as 
compared with the potential terms As in the present case, in equilibrium a first regime exists which is correctly 
described by naive mean-field theory. In the second one the critical temperature is given correctly by mean-field 
theory, but the critical exponents are not, and in the third one fluctuations dominate and mean-field theory has 
nothing to say. 

An important difference exists, however, between our nonequilibrium self-consistent solution for negative KPZ 
nonlinearities and equilibrium wetting in that the predictions for the former are for d = oo, while the three-regime 
behavior for the latter occurs only below the upper critical dimension d c = 2. In higher dimensions equilibrium 
wetting is known to be controlled by a Gaussian fixed point with trivial associated scaling [l[ . 

Higher-order moments also display interesting behavior. All moments mk starting from 2D /a 2 + 1 scale with the 
same exponent, while the usual scaling ~ m k obtains for k < 2D /a 2 + 1. This same behavior was observed in a 
mean- field study of nonequilibrium, MN1 complete wetting as reported in [Til ], and seems to be a common feature of 
multiplicative-noise controlled transitions. 

Additionally, it was found that the mean separation at the transition diverges as (h) ~ — ln(— b + b c ), with b c = 
for p < 2D /a 2 and b c < for p > 2D /a 2 . This implies that there is a single scaling regime in terms of h rather than 
three. That the behavior of the surface order-parameter is richer than that of the mean separation is seemingly a 
common characteristic of these systems. 

For positive KPZ non-linearities has also been investigated and the results show a single regime (n^ 1 ) ~ (— b) 1 ^, 
with higher moments scaling as m% ~ 77i fe . The mean interfacial separation grows logarithmically, (h) ~ mm ~ 
-ln(-6). 
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Hence, positive KPZ nonlinearities generate a trivial mean-field (or high-dimensional) behavior for nonequilibrium 
critical wetting, compatible with standard Gaussian scaling, which is analogous to the mean-field behavior of equi- 
librium critical wetting. On the contrary, negative KPZ nonlinearities behave in a rather intricate way, with highly 
nontrivial scaling including three different regimes even in mean-field (or high dimensions) approximation. 

In future work we will study how the rich phenomenology reported in this paper is affected by fluctuations, i.e. 
going beyond the mean field approximation. It would also be nice to have nonequilibrium critical wetting experiments 
to see whether our predictions can be observed in real systems. 
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FIGURE CAPTIONS 
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Fig. 1 (color online) Estimates of as a function of from power-law fits of (nop) vs b — b c after a numerical 
integration of Eq. ([5]) for MN1. Three regimes are found: the solid (blue) line corresponds to p < D/a 2 , the 
dashed (red) one to d/a 2 < p < 2D /a 2 , and the dash-dotted (green) one to p > 2D /a 2 . For ratios 2D /a 2 < p, 
b c is obtained from the (numerically computed) zeros of the corresponding parabolic-cylinder function. Close to 
the singularity p — D/a 2 the integration is cumbersome and the resulting errors large. Far from it, error bars 
are smaller than the symbols in some cases. 

Fig. 2 (color online) Estimates of /3 as a function of p from power-law fits of (nop) vs b — b c after a numerical integration 
of Eq. (O for MN2. The error bars are smaller than the symbols. 



